Elementary Excitations in Dimerized and Frustrated Heisenberg Chains 
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We present a detailed numerical analysis of the low energy excitation spectrum of a frustrated 
and dimerized spin S = 1/2 Heisenberg chain. In particular, we show that in the commensurate 
spin-Peierls phase the ratio of the singlet and triplet excitation gap is a universal function which 
depends on the frustration parameter only. We identify the conditions for which a second elemen- 
tary triplet branch in the excitation spectrum splits from the continuum. We compare our results 
with predictions from the continuum limit field theory . We discuss the relevance of our data in 
connection with recent experiments on CuGeOz, NaVzOs, and (VO)2P20i . 
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I. INTRODUCTION 

Low dimensional quantum spin systems have attracted 
considerable attention of theorists over the decades. 
Most of the interesting and fascinating features observed 
in these systems are pure quantum effects uniquely due 
to their low dimensionality. Peculiar properties of one- 
dimensional quantum antiferromagnets like e.g. exotic 
ground states or unconventional excitation spectra are 
not accessible to traditional methods like spin-wave or 
perturbation theory, but require the use of numerical or 
field-theoretical approaches. These methods are comple- 
mentary to each other and together with exact Bethe 
ansatz solutions of particular models they allow for a 
complete description of low dimensional quantum spin 
systems. In particular, the field-theoretical methods have 
been used successfully to predict the scaling behaviour of 
the one dimensional spin 5 = 1/2 Heisenberg model [nl, 
the existence of gapless and gapped phases in the 5 = 1/2 
frustrated Heisenberg chain and the existence of an 
excitation gap in the spin 5 = 1 Heisenberg chain ||. 
On the other hand only numerical methods allow to de- 
termine the critical value of frustration beyond which the 
gapped phase appears UH and to determine details of 
the ground state properties or the behaviour of the 
excitation gap itself 

Recently the interest in one-dimensional spin sys- 
tems has been particularly boosted by the discovery of 
various non-organic quasi one-dimensional compounds, 
in particular, the spin-Peierls materials CuGeO^ || 
and NaV 2 5 and spin ladder compounds like 

SrCu 2 3 , Sr 2 Cu 3 5 or possibly (VO) 2 P 2 7 [0-0. 
A common feature of these compounds is an excitation 
spectrum which is dramatically different from the spin 
5 = 1/2 Heisenberg chain. A remarkable fact about the 
Heisenberg chain is that its excitation spectrum consists 
of spin- 1/2 particles (spinons) [Q. Physically such exci- 
tations can be created only in pairs because upon flipping 



one spin the total spin projection is changed by AS Z = 1. 
Thus, in the Heisenberg chain the conventional magnons 
carrying spin 1 are deconfined into spin-1/2 spinons. In 
dimerized spin-Peierls compounds the excitation spec- 
trum is always gapped and the low lying excitations are 
triplets. In addition, a massive singlet branch may ex- 
ist above the triplet excitation branch in frustrated sys- 
tems. As it will be shown in the following, even a second 
triplet branch can appear below the continuum. Hence 
in these systems spinons are confined back into triplet 
magnons. The interaction between magnons can lead 
to massive singlet and triplet excitations below the con- 
tinuum [pi)| ]2Cf| and even a sequence of further massive 
excitations ||21|]. Furthermore the dimerized frustrated 
Heisenberg model can also describe the 2-leg ladder with 
frustration (22). 

In recent years the field-theoretical continuum-limit 
approaches were successfully used to st udy sp in-Peierls 
compounds and spin-ladder systems These 
studies show mechanisms for spinon confinement from 
the alternation of exchange couplings in spin-Pcicrls com- 
pounds and from the inter-chain coupling in spin-ladder 
systems. Although universal features of the physical 
system are usually properly captured in field theoreti- 
cal studies, important details governed by the physics at 
short length scales remain out of range for the applicabil- 
ity of these methods. Moreover, due to the perturbative 
nature of the continuum-limit approach, its predictions 
are less accurate in the physically more realistic strong 
coupling limit where details of the short distance physics 
are very important. Therefore, there is still a number of 
open questions motivating further theoretical studies of 
spin-Pcicrls and spin-ladder systems - especially in the 
strong coupling limit - by using exact methods. 

In this paper we present specifically a detailed numer- 
ical analysis of the low energy excitation spectrum of the 
5=1/2 antiferromagnetic Heisenberg chain with frus- 
tration and dimerization, as proposed in particular to 
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describe the magnetic properties of CuGeOs 
Hamiltonian reads, 



The 



where i denotes the sites of a chain with length L and 
Si are 5 = 1/2 spin operators. J > is the nearest- 
neighbor exchange coupling, a the frustration parameter 
from next-nearest neighbor coupling and S is the dimer- 
ization parameter. 

Besides its relevance to real spin-Peierls compounds 
the model is interesting purely from a theoretical point 
of view as far as it contains two independent mechanisms 
for spin-gap formation. At S = the model is charac- 
terized by a critical value of frustration a c which was 
accurately determined by numerical studies: a c = 0.2412 
|,§. For a < a c the frustration is irrelevant, the system 
renormalizes to the Heisenberg fixed point: the ground 
state corresponds to a spin liquid and the elementary 
excitations are massless spinons. At a — a c there is a 
transition into a spontaneously dimerized ground state. 
The spectrum acquires a gap and the elementary excita- 
tion is a massive spinon [|2j. On the other hand at any 
5^0 the singlet groundstate of the model is also dimer- 
ized with a gap in the spin excitation spectrum, but the 
elementary excitation is a magnon p|,p4[. 

For the special case 2a + 6 = 1 the groundstate of 
the spin Hamiltonian Eq. ([!]) is known exactly to be a 
product wavefunction of nearest neighbour singlet pairs 
|^5| , ^9| . This line in the (a, 5) parameter plane separates 
two distinct regimes: for 2a + S — 1 < the dominant 
peak in the static magnetic structure factor is at q* = 7r, 
while in the other case ir/2 < q* < it. In this latter 
incommensurate phase, q* continuously decreases from tt 
with increasing a and 5, and asymptotically approaches 
tt/2 0!. 

In a recent work pq] , the existence of a massive singlet 
excitation has been confirmed numerically for the Hamil- 
tonian Eq. (g) with a = 0.35 and 5 = 0.012. In addi- 
tion to the elementary triplet and singlet excitations and 
depending on the set of parameters (a,5) even another 
triplet excitation was found to split from the continuum 

IS 

In this paper we will analyze in detail the intriguing 
structure of the excitation spectrum for different frus- 
tration and dimerization parameters, and we compare 
our numerical data with the available results from field 
theoretical methods. The paper is organized as follows: 
In chapter II we summarize the procedure and the re- 
sults of the continuum limit field theory and outline the 
open questions inaccessible by these analytical methods. 
In chapter III we focus our attention on the singlet to 
triplet energy gap ratio, R(a, 8) — A s / At and show that 
R(a, S) only depends on a when 2a + 5 < 1; the field the- 
ory prediction R = is precisely realized for all i5 only 
when a = a c . In chapter IV we discuss the conditions 
for which a second triplet excitation branch may exist 



below the continuum. Finally, in chapter V we connect 
our results to recent experimental data on different spin 
chain compounds. 



II. THE CONTINUUM-LIMIT FIELD-THEORY 
APPROACH 

First insight into the structure of the excitation spec- 
trum of the model Eq. (Q) is obtained from bosoniza- 
tion and the continuum limit renormalization group ap- 
proach pjjj] . In terms of the continuum field cf>(x) the 
bosonized version of the initial spin model is the double 
sine-Gordon (SG) model 



Hbos = I d: 
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d4> 
dx 



M s cos(/?50) + M a coa(p a (/)) (2) 



where u is the spin wave velocity, (3 a — v87r, (3s = V 2n, 
M a = J(a — a c ), and Ms = S . 

is determined 



K = 



1 — (a— a c ) 



l+(a-a c ) ' 

The value of the critical frustration a c 
by the behaviour of the system at short distances i.e. 
it depends on nonuniversal parameters of the continuum 
limit theory. Therefore, differently constructed contin- 
uum limit theories give rather different values of this pa- 
rameter P, p9||3^ ] . The exact value of a c was determined 
only within numerical studies [^|,^| , but the non-universal 
parameters of the continuum limit Hamiltonian could be 
always chosen in such a way to ensure the proper value 
of the critical frustration a c . 

Contrary to the standard SG model (Hamiltonian || 
with only 1 'cosine term': Afcos(/30)), which is exactly 
solvable and well understood |Il) , the theory of the quan- 
tum double SG model is much less developed. However, 
in two limiting cases the model Eq. (||) reduces to the 
SG theory and provides exact knowledge about the char- 
acteristic properties of the system. 

Let us first consider the case 6 = (Ms = 0) corre- 
sponding to the frustrated Heisenberg chain. The be- 
haviour of this model is determined by the marginal in- 
teraction which is controlled by the frustration:/? = (3 a = 
^/8tt. For a < a c the interaction is irrelevant and the 
system scales to the Gaussian fixed point: elementary 
excitations are massless spinons. For a > a c the interac- 
tion is marginally relevant and the effective interaction 
renormalizes to large values. An exponentially small gap 
M* cx exp(— c/(a — a c )) is dynamically generated in the 
excitation spectrum, the field <p is ordered leading to a 
spontaneously dimerized ground state with the finite or- 
der parameter 



(Od) = -(^(-FrtSi-a -Si-Sf S i+ i)) 



(3) 



Thus, the elementary excitations in the massive phase of 
the frustrated Heisenberg chain are described by solitons 
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('kinks') of the quantum SG model with j3 = V8tt 
There are no soliton-antisoliton bound states in this case 
and the system is characterized by the only one scale - 
spin gap A = 2M*. Excitations above the given vacuum 
are created by breaking singlet bonds. Each broken bond 
gives rise to a pair of decoupled spins 1/2 on neighbouring 
sites. Once created, these isolated spins can propagate 
coherently along different sublattices and constitute ele- 
mentary excitations of the massive spinon type. 

We now consider the case 8 ^ and a — a c (M a = 
and Mg ^ ) The excitation spectrum of the SG model 
at /3 = (3g = V27T is exactly known |n| and at this point 
consists of soliton and antisoliton excitations with masses 
M s = M s = M and two bound states (breathers) with 
masses M\ — M and M% = \/ZM. The soliton excitation 
carries spin S z = 1, the antisoliton excitation S z = — 1, 
and the two breathers with opposite parity S z = 0. The 
lower energy breather mode is degenerate with the kink 
and anti-kink and these three excitations correspond to 
a triplet excitation branch in the original spin model lan- 
guage. The second bound state, in fact, has its counter- 
part in a spin singlet excitation || . These two modes are 
the only elementary excitations in this case and the ratio 
of their excitation gaps is exactly v3- 

The standard Renormalization Group (RG) approach 
gilll) to the double SG model Eq. @ is based on 
the fact that the critical dimensions of the two cosine 
terms arising from the smooth and staggered part of the 
exchanges, respectively, are different: 

dim cos(\f2~n<f)) = 1 , dim cos(v / 87T0) — 2 . (4) 

Thus, the 8 cos(-\/27r0) term is strongly relevant, while 
the J(a—a c ) cos(\/87r0) term is marginal. Therefore, the 
essential physics as determined by the relevant term is - 
at least for 8 -C 1 - similar to that of the above discussed 
SG model with /3 = \f2~n, and the marginal interaction 
leads to logarithmic corrections only. Therefore one as- 
sumes that the excitation spectrum of the spin-Peierls 
state consists of two excitation branches with gaps A t 
(triplet excitation) and A s = RA t (singlet excitation) 
with R slightly different from V3 due to the logarithmic 
corrections |24|,|3^]. However, since frustration and dimer- 
ization provide two principally different mechanisms for 
spin gap formation, interference between these interac- 
tions is non-trivial especially in the limit of strong initial 
interactions. A very sensitive tool to study these partic- 
ular effects is to explore the detailed structure of the ex- 
citation spectrum. Moreover, the exact excitation spec- 
trum of the dimerized and frustrated Heisenberg chain 
Eq. ([!]) provides another way for the determination of 
the critical parameter a c . Only for a = a c is the struc- 
ture of the excitation spectrum exactly the same as that 
of the SG model with (3 — ^/2~k. Therefore, a c is deter- 
mined from the condition R(S, a c ) = 

Due to the different critical dimensions of the cosine 
terms in the double SG model one may attempt, in a 
first approximation, to neglect the cos(-s/87r</>) term and 



to consider the usual SG model with the cos(v / 27r</()) term 
only. However, as we demonstrate in the subsequent 
chapters, we have to conclude from our numerical results 
that this commonly accepted procedure is not valid and 
the structure of the excitations is quite different. E.g. 
we find R = 2 in the absence of frustration a — which 
means that there is no long wavelength singlet excita- 
tion branch and the singlet excitation energy coincides 
with the edge of the continuum. Furthermore we obtain 
that for a < a c and small 8 the ratio R — A s /A t is a 
'universal' function of a alone. So although the contin- 
uum limit Hamiltonian is a proper description of the spin 
lattice model Eq. ([!]) the commonly adopted field theo- 
retical tools for the double SG model are not sufficient 
for a complete understanding of the excitation spectrum 
as we will show from our exact diagonalization data. 

III. GAP RATIO 

Our numerical study is performed using exact diago- 
nalization techniques with periodic boundary conditions 
for chains with up to L = 26 sites. As previously re- 
ported in Ref. the triplet and singlet gaps A t and 
A s have different finite size scaling behaviour. While 
At is a monotonically decreasing function of 1/L, A s is 
non-monotonic in 1/L and develops a minimum for a 
particular chain length which varies with the model pa- 
rameters. Thus, in order to extrapolate the values for A t 
and A s to the infinite chain limit, we need two different 
finite size scaling fit functions. For the triplet gap we 
have used the 3-parameter ansatz pq| , 

A t (£)=A t + ^exp(~) . (5) 

On the other hand, in order to account for the non- 
monotonic behaviour of A s we have chosen the 4- 
paramcter ansatz 

A s (L) = A s +exp('-AV^ + J B) (6) 

with A > and B < 0. Note that both fit functions 
proved to give an excellent agreement with density matrix 
renormalization group (DMRG) data with L of order 100, 
in particular in the region where the spin-spin correlation 
length is shorter than the chain length |2(||33| . 

In Fig.|l| we show the ratio R(a, 6) as a function of a 
for different values of the dimerization parameter (0.02 < 
S < 0.4) using the extrapolated values of A t and A s . 
Since we find that the width of the singlet dispersion 
is smaller than the width of the triplet dispersion, the 
necessary condition for the singlet excitation branch to 
split from the continuum over the whole Brillouin zone 
is R(a, 8) < 2. Fig.|l| clearly shows that the ratio de- 
pends on the frustration parameter in an essential way. 
For a < a c we obtain R(a, 8) > \/3 implying that the 
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cos(v / 87T</>) term in the double SG model is indeed rele- 
vant for all q/ci c and can not be omitted. We recall the 
field theoretical expectation that the deviations from V3 
should be logarithmically small. Furthermore, the exact 
diagonalization data show that the ratio R is for a < a c 
insensitive to the dimerization parameter S (to be more 
precise, for S > 0.05 within the accuracy of our finite 
size studies) which in the continuum limit field theory 
controls the relevant interaction in the double SG model. 
R = \/3 is indeed obtained at criticality a — a c for any 
finite <5, weak or strong. In the absence of frustration, 
a = 0, and for any dimerization we find R(a = 0, 8) = 2. 
Therefore, there is no well defined singlet excitation at 
q = (or equivalently q = ir) since A s = 2A t coincides 
with the lower edge of the continuum. 

In addition we find that (i) R(a, 8) is a universal func- 
tion of a for 2a + 8 < 1 and (ii) R(a, 8) depends on both 
parameters in the incommensurate phase for 2a + 8 > 1. 
We note that the deviations observed for weak values of 
8 (i.e. 5 = 0.05) in Fig.[l] result from a lack of preci- 
sion in the determination of the extrapolated A s . High 
precision is lost when the spin-spin correlation length be- 
comes comparable to or longer than the chain length for 
small 8. On the other hand the comparison with DMRG 
calculations for L < 100 shows that the ansatz for the 
scaling of the singlet gap A s , Eq. (^), remains very accu- 
rate even in the region of long correlation lengths (small 
gap). To demonstrate the consistency with the diagonal- 
ization data we have added in Fig.|l| DMRG data points 
calculated for 8 = 0.02 [f|. The DMRG data support 
the observation that the ratio R(a, 8) depends on a only 
in the commensurate phase, i.e. as long as 2a + 8 < 1. 
We have indicated in Fig.|l| (with a cross) the data points 
for which the parameter pairs (a, 8) belong to the in- 
commensurate phase. We emphasize that for these sets 
of parameters the finite size effects are extremely small 
and the values of i?(a, 5) are thus very accurate. 

A similar behaviour has recently been found by Yoko- 
yama et al. for the leading 8 power law dependence of 
the triplet excitation gap A t oc 8 1 , where the exponent 

7 is a monotonic continuous function of a |27|| . In this 
work it was shown that the Cross-Fisher value 7 = § p4| ] 
is realized only for a = a c , but 7 ^ I for a < a c . 

IV. A SECOND TRIPLET BRANCH 

Another peculiar observation is made when the param- 
eters (a, 8) are increased towards stronger dimerization 
in that a second triplet excitation branch splits from the 
continuum. In analogy to the singlet excitation this sec- 
ond triplet may be interpreted as a bound-state between 
a triplet and a singlet excitation. In order to investi- 
gate this new feature we fix the dimerization parameter 

8 = 0.2 and discuss the low energy spectrum as a function 
of frustration in the subsector of total momentum q = 0. 
For this purpose we have show in Fig.Sa the energies of 



the 3 lowest excited states, extrapolated to L — > 00, as a 
function of a. In this figure the lower edge of the contin- 
uum at 2A t is indicated by a continuous bold line. We 
observe that for a < 0.28, there are only two well defined 
excitations below the continuum, one triplet and one sin- 
glet as discussed above. The energy of the next excited 
state is found to scale to the lower edge of the contin- 
uum with L — > 00. In contrast to the triplet, the singlet 
excitation gap remains almost constant with changing a. 
However, we observe that for a > 0.28 another triplet ex- 
citation T 2 splits from the continuum (i.e. E? 2 < 2A t ). 
This is in agreement with the results reported previously 
in Ref. 

In order to verify that the triplet T2 is indeed a well de- 
fined elementary excitation at q = 7r only when a > 0.28 
(for 8 — 0.2), we have also evaluated its spectral weight 
W = I (2^2 1 1 0) 1 2 versus 1/L in the dynamical structure 
factor S(tt, lu) for different frustration parameters a (see 
Fig.^b) where |0) is the groundstate wave function and 
6 = S 2 (tt) = (l/-y/Z) £ z exp(i7r Z)Sf. We observe that 
for a < 0.28, W has a strong size dependence, and the 
data indicate that W — > with L — > 00. However, for 
a > 0.28 the curvature of W(l/L) changes indicating 
that the weight W(l/L) scales to a finite value consistent 
with the identification of T2 as an elementary excitation 
below the continuum. 

By evaluating the excitation spectrum in different mo- 
mentum sectors we find that the triplet T2 splits from 
the continuum first at q = J. In Fig.||a, we have plotted 
the ratio, 

g(S = l,f,2)-£(S = l,f,l) 
Ri = Xt (7) 

versus a where E(S, q, n) is the n-th energy level in the 
subsector with total momentum q and spin S extrapo- 
lated to L — > 00. Ri < 1 implies that the excitation is 
split from the continuum at q = ^. FigJ^a shows that 
for large enough dimerization (e.g. i5 = 0.2) the triplet 
T2 is well defined at q — ^ for any a. However, as we 
reduce the dimerization parameter (e.g. to S = 0.1), we 
observe that stronger frustration is needed in order to 
separate T2 from the continuum. In fact, for 8 — 0.1 and 
a < 0.12, Ri « 1, i.e. the energy of the second triplet ex- 
citation scales to the lower edge of the continuum. Thus 
our data suggest that in the absence of frustration, there 
is a minimal value 8 m i n > 0.1 above which the triplet 
T2 appears. The fact that i?i ~ 0.98 instead of 1, for 
8 = 0.1 and a < 0.12 shows that this excitation is not 
well defined, the finite size analysis is not valid for excita- 
tions in the continuum. The kink in R± clearly indicates 
the minimum value of the frustration for which the ex- 
citation separates from the lower edge of the continuum. 
Similarly, the value a = 0.28 (see Fig. |2^,) is the min- 
imal frustration a m i n (8 = 0.2) for the second triplet to 
be a well defined excitation branch over the entire Bril- 
louin zone. In other words, for a given 8, the momentum 
range for which the triplet T2 is split from the continuum 
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is centered around q = tt/2 increasing continuously with 
increasing a. As an example, we show in Fig.^the disper- 
sion of the 3 lowest excitations for a = 0.2 and 8 = 0.2, 
for which the triplet X2 is split from the continuum in 
the momentum range Air/5 > q > q m in ~ 7r/5. 

In order to determine the (5-dependence of a m i n , we 
have plotted in Fig.||b the ratio 

R 2 = E{S = 1,«,2)-Eq (8) 
At 

versus a for different values of 8. The excitation T2 is 
split from the continuum when R 2 < 2. Fig.^b shows 
that a m in increases when 8 is reduced, for instance 
a mm (0.2) « 0.28 and a mi „(0.05) « 0.32. We have plot- 
ted in Fig.|]c a m i n as a function of 8 for 8 < 1. Our results 
suggest that a m i n is within a good approximation a sim- 
ple linear function of the parameter 8, a m %n ~ (1 — 8)/3. 
Thus, in the special limiting case 8 — ► 1 , and restrict- 
ing ourself to the commensurate region we observe that 
ct-min —> 0. Note that the result can not be extended at 
the special point 8=1, since the spectrum at this point 
consists of a discrete set of eigenvalues E n = n * At, thus 
there is no continuum in the vicinity of 2 * A t . How- 
ever, for 0.5 < 8 < 1, the spectrum consist of mini-bands 
with finite band-gaps, thus the notion of continuum has 
a meaning in the vicinity of 2 * At . Furthermore our re- 
sult suggests that the spectrum is completely dense (no 
band gap), when 8 < 0.5. The other interesting point 
corresponds to the limit 8 — > 0. Indeed, if the linear 
approximation is still valid in this limit, our data shows 
that a mvn —►1/3. 

V. DISCUSSION OF EXPERIMENTAL DATA. 

We now try to connect these results to some recent ex- 
perimental data on CuGeOs, Nal^Os and (^0)2^207. 
CuGe0 3 § and NaV 2 5 are spin-Peierls sys- 

tems with transition temperatures T$p ~ UK and 35-ftT, 
respectively. Previously it has been shown that the 
magnetic susceptibility of CuGeO^ in the uniform phase 
can be accurately reproduced by a frustrated Heisenberg 
chain model with j37|]; from this fit the frustration was 
estimated to be q ~ 0.35, close to the previous estimate 
of Riera and Dobry j38|. With an exchange coupling 
J sa 160K and a — 0.36 the experimental value for 
the gap A t = 2.1meV as determined by inelastic neu- 
tron scattering |}6| is obtained within the frustrated and 
dimerized Heisenberg chain model with a dimerization 
8 = 0.012 |26| ], With this parameter set a singlet-triplet 
gap ratio i?(0.35, 0.012) « 1.50 follows (see FigJ|). How- 
ever, the experimental ratio is R eX p = 1-72 w \/3; ac- 
cording to Fig.|l|, this would corresponds to a w a c , a 
value for a which was previously proposed by Castilla et 
al. §. 

There arc different possibilities how to resolve this 
quantitative problem. First, the effects of interchain cou- 



pling along the crystal 6-direction in CuGeOz are exper- 
imentally well established but little is known from theo- 
retical studies |H . Second, the dynamics resulting from 
the spin-lattice coupling is expected to play an important 
role. Indeed, including dynamical phonons will renormal- 
ize the spin excitation spectrum and as a direct conse- 
quence the ratio R ]4C| ]. 

Recently, it has been proposed that the low temper- 
ature phase of the compound iVaV^Os is well described 
by an unfrustrated (a ~ 0) dimerized chain model [ jTo| . 
Following the same procedure as used for CuGeOs, the 
dimerization was estimated to be of order 8 « 0.048 JIT| . 
If we assume that this is indeed the case then - accord- 
ing to our results - we expect a singlet-triplet gap ratio 
R(0, 0.048) = 2, i.e. there is no well defined long wave- 
length singlet excitation below the continuum. As a con- 
sequence, we expect no magnetic Raman response at fre- 
quencies below the lower edge of the continuum at 2A t . 
However, in a recent Raman scattering experiment on 
NaV20s Q it has been observed that the energy of the 
lowest excitation in the dimerized phase is at 64cm -1 ; for 
this compound the measured triplet excitation gap from 
inelastic neutron scattering is A t « 59cm -1 [^lj. Thus 
the energy of the lowest excitation in the Raman spectra 
is very close to the singlet-triplet gap A t . Furthermore, 
it was observed that this feature remains unchanged for 
different photon polarization geometries suggesting that 
this excitation is probably not of magnetic origin. An 
alternative origin is a charge excitation due to a local 
charge transfer from a V A+ to a V 5+ chain. All these 
observations indicate some important qualitative differ- 
ences between the spin chain compounds CuGeO% and 
NaV 2 O b . 

Finally we discuss the case of (^0)2^207 ~ a spin 
chain compound whose magnetic properties have given 
rise to controversial interpretations. This compound has 
been initially considered as an ideal realization of a two- 
leg antiferromagnetic Heisenberg spin ladder . Subse- 
quently, however, it has been shown that the susceptibil- 
ity data on polycristalline material could be well fitted 
by either a ladder or an alternating chain model p4[ . 
Early neutron scattering data indicated a spin gap of 
about 50^ and supported the two-leg ladder picture jl5) . 
Recently, inelastic neutron scattering experiments have 
been performed on a collection of many oriented single 
crystals jL5|. The observed features like e.g. the strong 
dispersion in the rung direction lead to the conclusion 
that the data are not consistent with a spin ladder de- 
scription and that (V 0) 2 P20'j is better described as an 
alternating spin chain directed in the rung direction. Re- 
markably, in addition to the lowest triplet another well 
defined higher energy triplet excitation below the contin- 
uum was observed near q = tt. If we assume that this 
compound is well described by an unfrustrated (a = 0) 
dimerized chain with a dimerization of order 8 ~ 0.1 (i.e. 
an alternating exchange coupling in the rung direction), 
then we find a contradiction with our present results. In 
fact, for 8 ~ 0.1 and a = 0, there is no well defined 
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second triplet at q — n. The microscopic structure of 
{VO^P^Ot, however, does not allow to identify any ob- 
vious superexchange path which can lead to frustration. 
At this point we can only point out that both pictures - a 
two-leg ladder or an unfrustrated alternating chain - do 
not properly describe the low energy excitation spectrum 
in this compound. The correct modeling of the magnetic 
properties of (yO^F^C^ remains an unresolved problem. 



VI. SUMMARY AND CONCLUSION 

In this paper we have studied the elementary excita- 
tions of frustrated and dimerized Heisenberg chains. We 
have shown that in the commensurate dimerized phase 
the singlet-triplet gap ratio R(a, S) only depends on the 
frustration parameter for arbitrary values of the dimer- 
ization (5 > 0.02). The magnitude of this ratio gives a di- 
rect measure for the strength of the frustrating exchange 
coupling. Without frustration R(a = 0, 6) = 2 implying 
the absence of a long wavelength singlet excitation below 
the continuum. We have shown that away from critical- 
ity a a c the frustration term plays an important role 
for the low energy excitation spectrum. Therefore, the 
commonly adopted procedure of disregarding the frus- 
tration term in the bosonized continuum limit Hamilto- 
nian for a < a c is at least questionable since it is not 
marginally irrelevant. However, for a = a c , our exact di- 
agonalization results give precisely R(a c ) — \/3 in perfect 
agreement with the field theoretical treatment of the SG 
model. Furthermore, we have demonstrated the condi- 
tions for the existence of another triplet excitation branch 
which we interpret as a bound state of elementary triplet 
and singlet excitations. Depending on the dimerization 
strength the second triplet excitation branch is observ- 
able at q = (or q = tt) only if the frustration parameter 
is large enough. As we have discussed we consider our 
results relevant for the magnetic excitations in the spin 
chain compounds CuGeOz, NC1V2O5, and {yO)iP20i. 
For the latter material we conclude from our results that 
neither a two-leg ladder nor an alternating, unfrustrated 
Heisenberg chain model captures the experimental facts 
on the excitation spectrum in (V 0)2^207 . 
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FIG. 3. (a) Ri and (b) R2 as a function of a calculated 
for different values of the parameter 8. For the definitions of 
the ratios Ri and R2 see Eqs. (Q) and (^), respectively. In fig 
(c), we have plotted a m in as a function of 8, the dashed line 
is a linear fit of the data. 



FIG. 4. Dispersion uii(q) = Ei(q) — Eo of the three lowest 
excitations for a — 0.2 and 8 = 0.2. The full symbols corre- 
spond to the triplets T\ and T2 and the open symbols to the 
singlet Si. The continuum is indicated by the shaded area. 



FIG. 1. Extrapolated singlet-triplet gap ratio R(a, 8) — 
A s /At vs. frustration a for different dimerization parame- 
ters 8. The symbols marked with a cross have been calcu- 
lated in the incommensurate phase. ED indicates exact di- 
agonalization results extrapolated to the infinite chain limit. 
In addition we have added data points from density matrix 
renormalization group (DMRG) calculations J33I. 



FIG. 2. (a) Energy of the 3 lowest excitations vs. frustra- 
tion calculated at fixed dimerization 8 — 0.2 and momentum 
q = (resp. n). The continuous line indicates the edge of 
the continuum at 2 At. Ti,2 are triplet excitations and Si is 
the lowest singlet excitation, (b) Spectral weight W of the 
triplet T2 as a function of 1/L (L is the size of the system) for 
8 = 0.2 and different values of a. The full symbols correspond 
to cases for which W scales to finite values in the infinite chain 
limit. 
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Singlet-Triplet gap ratio vs. frustration (Q=0). 
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0.080 

Weight 

0.070 
0.060 
0.050 
0.040 



G- 


^Da=0 


□- 


-B a=0.1 


0- 


-0 a=0.2 


A 


^4 a=0.3 




a=0.4 



8=0.2 




0.030 



0.020 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 ' 

0.000 0.010 0.020 0.030 0.040 0.050 0.060 0.070 0.080 

1/L 



R, 



1.08 
1.04 
1.00 
0.96 
0.92 
0.88 
0.84 
0.80 











A A 5=0.2 






■ ■ 5=0.1 












0.76 
0.00 



0.10 



0.20 



0.30 



a 



Triplet-Triplet ratio as a function of alpha. 




1.40 ^ ' ^ ' ^ ' ^ ' 

0.10 0.20 0.30 _ 0.40 0.50 



0.40 



a 



mm 



0.30 - 



0.20 - 



0.10 - 



0.00 




